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Abstract 

In the first part of tliis introductory review we outline tfie develop- 
ments in pliotonic band gap materials from tlie pliysics of pliotonic band 
gap formation to the fabrication and potential applications of photonic 
crystals. We briefly describe the analogies between electron and photon 
localization, present a simple model of a band structure calculation and 
describe some of the techniques used for fabricating photonic crystals. 
Also some applications in the field of photonics and optical circuitry are 
briefly presented. In the second part, we discuss the consequences for the 
interaction between an atom and the light field when the former is em- 
bedded in photonic crystals of a specific type, exhibiting a specific form 
of a gap in the density of states. We first briefiy review the standard 
treatment (Weisskopf-Wigner theory) in describing the dynamics of spon- 
taneous emission in free space from first principles, and then proceed by 
explaining the alterations needed to properly treat the case of a two-level 
atom embedded in a photonic band gap material. 



1 Introduction 

At the microscopic level, ordinary matter exhibits behaviour analogous to light 
waves. When wave-like electrons scatter off ions in crystalline materials, con- 
structive interference between different trajectories can cause electric currents 
to flow. Conversely, disorder in such crystals can hinder electrical conductivity, 
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and for some energies the electrons become localized in space, thus preventing 
their free flow in the form of electrical currents . Although such ideas have been 
known since the 1960s, it is only recently that physicists have begun to ask 
whether similar effects can result in the localization of light in a corresponding 
"photonic crystal" . 

The purpose of this article is twofold. In the first part we outline the devel- 
opments in photonic band gap (PBG) materials from the physics of photonic 
band gap formation to the fabrication and potential applications of photonic 
crystals. In the latter, we specifically focus on the state of art 3-D structures 
exhibiting a full band gap in their electromagnetic field density of states, be- 
cause as we believe, these are of great technological and fundamental interest. 
In the second part, we discuss the consequences for the interaction between an 
atom and the light field when the former is embedded in dielectric materials, 
mostly photonic crystals, which exhibit gaps in the density of states. We first 
review the standard treatment to describe the dynamics of spontaneous emis- 
sion in free space and proceed by explaining the alterations needed to properly 
illustrate the dynamics of two-level systems embedded in a PBG material. 

2 Developments in photonic crystals 
2.1 Prom electrons to photons 

A normal crystal is a periodic array of atoms which scatters and modifies the 
energy momentum relation of electrons, whereas a photonic crystal is an ordered 
inhomogeneous medium characterized by a spatially periodic dielectric constant, 
with lattice parameter comparable to the wavelength of the light p. Fig. ^ 
shows one such microstructure in which a complete bandgap at 1.5 microns has 
been observed. 

Strictly speaking, such a structure has no allowed electro-magnetic (EM) 
modes in the forbidden range of frequencies or in other words the density of 
states of the propagating photon modes is zero (see Figs. El El- By contrast in 
free space (a cavity of infinite volume) the density of modes varies as and 
exhibits no gap. In a cavity of finite volume the density of states is substantially 
modified for frequencies close to the cavity cut off. Below the cut off the cavity 
sustains no modes at all, and near above the cutoff, the density of states can be 
increased to the continuum case. The vanishing of the density of propagating 
photon modes within a PBG means that, for the frequency range spanned by 
the gap, linear propagation of EM waves is forbidden in any direction in the 
PBG material. Thus, light incident on a PBG material with a frequency from 
the gap region will be backscattered from the material, independent of the angle 
of incidence. Strong suppression of transmission with an associated peak in the 
reflectivity at the characteristic frequencies is then an experimental signature of 
a photonic band gap 0. 

To give a further insight on the physics of light localization we will exploit the 
analogy between electron and light a bit further. As is well known an electron 



2 



Figure 1: A 3-d silicon photonic crystal exhibiting a complete band gap for 
incident light at 1.5 microns. It is constructed by growing silicon inside the 
voids of a opal template of close-packed silica spheres that are connected by 
the small "necks" formed during sintering, followed by removal of the silica 
template-see Sec. 2 (courtesy of S. John's group, University of Toronto). 




frequency frequency 

Figure 2: The density of modes Figure 3: The density of modes of the 

of the EM field in free space. EM field in a PBG structure with a 

broad gap. 
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in a disordered solid is described by the following Schrodinger equation: 



- ^ W(x) + V^(x)V'(x) = i?^(x) , (1) 

where m* is the electron's effective mass and V{'s) is a potential that varies 
randomly in space. For sufficiently negative energies E, the electrons may get 
trapped in regions where the random potential is very deep. The rate at which 
electrons tunnel out of the deep potentials depends on the probability of finding 
nearby potential fluctuations into which the trapped electron can tunnel. This 
rate increases as the electron energy increases. 

In the case of a monochromatic EM wave of frequency uj propagating in an 
inhomogeneous but non-dissipative dielectric medium, the classical wave equa- 
tion for the electric field is 



I ap^ 



where P is the polarization of the medium. Assuming now propagation in a 
linear medium, this leads to 

- V^f - --2-e/i„ct(x)£ = en-^f , (3) 

where S is the slowing varying field amplitude. Also the total dielectric constant 
has been separated into its average value cd and a spatially fluctuating part 
Sfiucti'^)- In a lossless material the dielectric constant e(x) is everywhere real 
and positive ^ and plays a role analogous to the random potential V^(x) in the 
Schrodinger equation. It scatters the EM wave. 

Comparing now the above two equations ^ and Q we observe the following 

2 

differences: Firstly, the quantity cdtt which plays the role of an energy eigen- 
value, is always positive, which precludes the possibility of elementary bound 
states of light in deep negative potential wells. 

Secondly the mode frequency lu multiplies the scattering potential e fiuct (x) and 
in contrast to an electronic system, where localization is increased by lower- 
ing the electron energy, lowering the photon energy leads to a complete dis- 
appearance of scattering. In addition, looking at the opposite high-frequency 
limit, geometric ray optics becomes valid and interference corrections to opti- 
cal transport becomes less and less effective. These simply mean that in both 
cases the normal modes of the EM field are extended, not localized. Finally, 
the condition that eo + Sfiuct > everywhere translates into the requirement 
that the energy eigenvalue will always be greater than the effective potential 
|^e/iMct(x)|. Therefore, unlike the familiar picture of electronic localization, 
what we are really seeking is an intermediate frequency window within the pos- 
itive energy continuum that lies at an energy higher than the highest of the 
potential barriers. For those frequencies the interference between the incident 

^AU the dielectric materials described in this article will assumed to be completely lossless. 
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and the scattered EM waves will be exactly destructive and this will allow no 
free propagating photon modes to exist in the structure. 

In the following section, we quantify these ideas in a more rigorous way and 
present the procedure followed to find the allowed modes as a function of the 
frequency, the dispersion relation^ in a simple 1-D periodic dielectric. 



2.2 Calculating the dispersion relation 

The dielectric constant of any photonic crystal can be expressed as 

e(x) = ei) + e/;„ct(x) . (4) 

In this case we will assume that our structure exhibits periodicity in one dimen- 
sion and is homogenous in the other two. More specifically 

+00 

efiuct{x) = ^ u{x-nL), (5) 

n— — 00 

with L being the lattice constant and 

u(x) ^ I \x\ <a .g, 

^ ' \ 0, otherwise ' ^ ' 

As discussed in the previous section, the propagation of a monochromatic EM 
field in an inhomogeneous, non dissipative dielectric medium is governed by the 
following equation 

2 2 
- V^f + V(V • £) - ^efiucti^)£ = ^D^£. (7) 

where £ is the field's ampitude. Setting = 1, the above equation reads 
(V • £ = 0, free charge is zero) 

2 2 

- V'^Eix) - '^efiuct{x)£{x) = ^^Six) . (8) 

Setting $(a;) = -^efiuct{x) we get 

^2 

- \/^£{x) + ^{x)8{x) = -^S{x) , (9) 

where the potential <I>(x) is basically a sequence of "potential barriers" of width 
2a (see Fig. 01. The objective is to solve Eq. ^ for the potential $(a;) 

Restricting ourselves to a unit cell of the crystal, the solution can be ex- 
pressed as follows 

j^^tkx ^ ^-tkx ^ x < -a 
£{x)={ Ce'^''' + De-'^'"", \x\ < a , (10) 
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Figure 4: The scattering potential for electrons in a solid (top figure) and for 
photons in a random dielectric medium (bottom) [2]. The effective scattering 
potential for photons is (w^/c^ )£/;„£*, where ejiuct is the spatially varying part 
of the dielectric. The electron (top) can have a negative energy, and can be 
trapped in deep potentials. In contrast, the eigenvalue {uj'^/c^)€d of the pho- 
ton (bottom) must be greater than the highest of the potential barriers if the 
dielectric constant is to be real and positive everywhere. 
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Figure 5: The periodic potential. L is the lattice constant 



where k = ui/c and k' = nuj/c. According to the Floquet theorem, £{x) should 
obey the following equation 



E{x + L) = e''''^E{x), (11) 



and the derivative 



d£{x + L) ^ ^,kL d£{x) ,^2) 
dx dx 

In addition the boundary conditions are that both £{x) and d£{x)/dx should 
be continuous at a; = ±a. Applying all these conditions, we calculate the 
expansion coefScients and also the following transcendental equation describing 
the relation between oj and k: 

ii,T\ ,'2riauj buj r? + 1 . 2nauj . hw 

cosikL) — cost jcosl — 1 sin[ )sin{ — ) , [16) 

c c 2n c c 

which for b — 2na can be inverted analytically providing the dispersion relation 
for the ID crystal, 

c AncoslkL) + [l - n"^) ^ 
Wfc = - — arccos[ ^— ^ — ^ '-] . (14) 

Ana (1 + 

This dispersion relation leads to gaps at fc = 2{n+i)a ^'^^ "-"^"^ integer values of 
m (see Fig. EJ. The lowest gap is centered at the frequency ujgap = 7rc/(4na) 
which for the case b = 2na equals tt/L (case shown Fig. 

In the more general case of a realistic 3-D photonic crystal, the usual ap- 
proach is using the Bloch-Floquet theorem to expand both the field amplitude 
and the dielectric constant in plane waves whose wavevectors are reciprocal 
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Figure 6: The dispersion relation for the 1-D isotropic model shown in Fig. |31 
As illustrated, gaps are formed at rrni jL for integer values of m. 




Figure 7: The band structure calculation of a realistic 3-D photonic crystal(the 
inverse opal structure of Fig. 1) according to theoretical calculations. The 
shaded region corresponds to the gap predicted (width 5.1%) and partially ob- 
served by the Toronto group 
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lattice vectors. 

^i(r)=EVlele'(''+''^ (15) 

G,l 

where e/ is the polarization vector and G are reciprocal lattice vectors. Using 
this expansion in Eq. Q, the problem is reduced to the solution of a sys- 
tem of linear equations on bg i the solution of which provides the allowed mode 
frequencies for a given crystal (Fig. (Tj). We should note here that for more 
complicated 3-D structures the solution of the corresponding eigenvalue equa- 
tion proves to be cumbersome and highly efficient computational techniques are 
usually required, see for example Refs. |31 1^ ITIIFII^ [TUl [TT| . 

2.3 Fabrication of photonic crystals 

In order to realize a 3-D structure with a full band gap for the propagation 
of light for some specific frequencies, we need not only to show that a spe- 
cific geometry could in principle exhibit a band gap but also that the specific 
micro-structure is amenable to microfabrication in the lab |12j. To illustrate 
the methods followed by various groups in this field we will start with the sim- 
plest case, the traditional multilayer film. The latter is the simplest dielectric 
structure where one can observe inhibition of the linear propagation of a EM 
wave. It is relatively easy to construct by assembling together dielectric layers 
with alternating high and low refractive indexes. 

To go further than that, i.e. to create a band gap in the propagation of an 
EM wave in two dimensions, we need something more sophisticated than simply 
contrast in the index of refraction. We need to find a specific geometry that will 
provide a full band gap for a range of frequencies in two dimensions. For this 
we need to keep in mind that the EM field consists of two type of modes. The 
transverse-electric (TE) and transverse magnetic (TM) ones. To achieve a com- 
plete band gap for all polarizations, the corresponding bands should not only 
exist but also overlap. TM band gaps are favored in a lattice of isolated high 
index regions as in an array of dielectric columns in air. In contrast, TE band 
gaps are favored in the inverse structure as in an array of air columns (veins) 
drilled in a dielectric substrate. Therefore to achieve a full band gap, we need 
somehow to reconcile these seemingly contradictory conditions ^.1 . A structure 
satisfying this, is a triangular lattice of low index columns (air) inside a high in- 
dex medium (silicon), see Fig. |S1 More specifically for the case where the radius 
of the columns is large enough, the spots between the columns behave like lo- 
calized regions of high refractive index material and thus the above requirement 
should be satisfied. It was predicted that for r/a ~ 0.48 for the ratio between 
the radius of the air columns to the lattice constant, a complete photonic band 
gap should form exhibiting a gap-midgap frequency ratio of about 19%. The ex- 
perimental verification was provided by Gruning et el. 13 , who constructed the 
crystal using an electrochemical technique to etch out columns in a silicon sub- 
strate. Their subsequent measurements verified a band gap at A = 5/im. Also 
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Figure 8: Macroporous silicon forming a 2-dimensional triangular lattice taken 
using scanning electron techniques ^S]- (a) A porous silicon bar of width 200 fj.m 
and height 75/iTO. (b) A tenfold magnification of the inset shown in (a), (c) A 
tenfold magnification of the inset shown in (b) . The lattice constant is of the 
macropore array is 2.3fini, the pore diameter 2.13/zm, and the thinnest parts of 
the pore walls 170nm. 

at the same year Krauss et al. constructed 2D polarization-sensitive photonic 
band-gaps with wavelengths in the range 800 - 900 nm . 

However, to achieve complete inhibition of light in all directions, we need 
a 3D structure which exhibits periodicity in all three directions. Out of the 
plethora of different geometries, the one that was found, in theory first [T5] 
and experimentally later [W\ to support a full 3-D band gap is the diamond 
lattice ^. The structure initially suggested by the Iowa group consists of ei- 
ther dielectric spheres in air or air spheres embedded in a dielectric medium. 
Yablonovitch managed to implement the latter by mechanically drilling cylin- 
drical holes through a dielectric block (with n ^ 3.6). The points where his 
tunnels coincided formed a diamond-hke structure (see Fig. EJ. In spite of the 
apparent simplicity and initial success of the method in demonstrating the exis- 
tence of a band gap in the microwave region, the application of the same method 
for optical waves proved to be of great difficulty. In this case, the spheres and 
consequently the drill had to be in the micrometer region! We remind here that 
in order for the TM and TE gaps to overlap, both constituent materials (air 
and dielectric) have to be topologically interconnected and also a large contrast 
(at least 3) between the two is essential. These constraints limited the use of 
microengineering fabrication techniques (electron beam and X-ray lithography 
[T7| which proved to be very successful in longer wavelengths ^ . 

^We note here that proposed common semiconductors such as siUcon and germanium also 
follow diamond geometry. 

•^We note here the ion drill proposed by the Yablonovitch group which led to the fabrication 
of structures with band gaps in the near-infrared (1.1-1.5/im). The problem in this case was 
that only a few unit cells could be reliably produced. 
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Figure 9: The first prototype structure predicted to exhibit a large and robust 
3D PBG structure designed by Yablonovitch jT^ . It consists of an overlapping 
array of air spheres arranged in a diamond lattice. It was created by drilling an 
array of criss-crossing cylindrical holes into a bulk dielectric of refractive index 
n = 3.6 where band gaps of the order of 20% where demonstrated. (Picture 
courtesy of the Toronto group) 

A different approach was the layer stacking technique proposed by Ho et 
al. 5H| (see Fig. and implemented by Lin et al. where a gap at 1.5 
microns was reported and the seven layer crystal exhibited 1% transmission 
(it was believed to drop to 0.1 % with ten layers). Along similar lines using 
a combination of electron beam lithography and reactive ion etching Noda et 
al. managed to stack semiconductor rods with micrometer dimensions |20| . 
They reported 99.9 % attenuation between 6 and 9 jim using eight layers. We 
note here that both methods are promising for any large scale technological 
application (see next section) as they are characterized by low cost and high 
reliability. ^ 

To overcome the difficulties of 3-D sub-micron engineering it has also been 
proposed to utilize systems that tend to self assemble themselves into various 
geometries. The most promising one proved to be colloidal crystals and ar- 
tificial opals. Colloidal particles have been synthesized by materials such as 
latex and Si02 in the range of a few nanometres to a few hundred micrometers. 
A suspension of colloidal microspheres, with a typical concentration of 10^° 

*The MIT group led by Joannopoulos has recently proposed to build a 3-D crystal by 
stacking 2-D ones (see Fig llOI . They predict the opening of a full 3-D band gap. 
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Figure 10: Computer reproduction of a novel 3-D photonic crystal suggested by 
the Joannopoulos group [21] ■ The structure could be fabricated layer by layer 
following a three layer period thus allowing the use of traditional lithographic 
techniques, along with a high degree of control in placing defects (waveguides, 
cavities, and other optical components) in the crystal. The layers consist of al- 
ternating stacks of the two characteristic types of 2-D (or slab) photonic crystals: 
dielectric rods in air and air holes in dielectric. 

particles/cm'^, can residue under gravity into a cubic-closed-packed structure 
with size of the order of 1cm. ^ 

In spite of the advantage of producing inherently 3-D structures com- 
pared to the 2-D ones using lithographic techniques, it appears to be extremely 
difficult to achieve the required refractive index contrast and interconnected- 
ness for a full 3-D band-gap to open. On the other hand, if they are used as 
templates to fabricate inverse opals (Fig^l, i.e., cubic-closed-packed lattices 
of air bubbles on a dielectric matrix (silicon or GaAs) 0] I23L , near visible 
photonic band gaps with gap to mid-gap frequency ratio of about 10% were 
predicted (see Fig. [7|). We note here that the void regions that are left behind 
after the etching of the original template (Swiss cheese structure) will allow 
the injection of atoms or dye molecules thus making quantum optical type of 
experiments possible.^ In spite of the initial success of the inverted opals in 
providing rather large band gaps in infrared region, problems of gap instability 
to disorder effects may prove to be difficult to overcome. A solution to this 
was recently proposed by S.John's group in University of Toronto [211 ■ They 
proposed an alternative photonic crystal architecture consisting of square spiral 
posts in a tetragonal lattice (Fig. [T^ . It seems that this structure could exhibit 
a quite large and robust 3D PBG occurring between the fourth and fifth EM 

^Their crystalline structure resembles that of natural opals. 

^In Ref. 1251 a similar type of experiment for Bell Inequality test is proposed. 
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Figure 11: The woodpile structure fabricated by repetitive deposition and etch- 
ing of multiple dielectric films, a, Sketch of the layer-by-layer 3-D photonic band 
gap structure. It is made from layers of one-dimensional rods with a stacking 
sequence that repeats itself every four layers, a unit cell, with a repeat distance 
of c. b. Plot of the computer calculated photonic density of states assuming 
that the refractive index of the rods is n = 3.60, the filling fraction of the 3-D 
structure f = w/d = 0.28, and c/d = 1.414. A complete photonic bandgap exists 
from 0.46 c/ a to 0.56 c/ a, where c is the speed of light in vacuum |17| . 
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bands and is also amenable to large-scale microfabrication using glancing an- 
gle deposition (GLAD) techniques. Spiral post lattices with microscale features 
have previously been synthesized using the GLAD method. In this technique, 
complex 3D structures can be fabricated by combining oblique vapor deposition 
and precisely controlled motion of a two-dimensionally patterned substrate. In 
addition, their square spiral posts can serve as templates for growing PBG mate- 
rials from an even larger range of materials. In this high refractive index 
material may be infiltrated to fill the void regions between the posts, with the 
posts subsequently removed by some selective etching process, leaving behind 
an "inverted structure." Their calculations have shown that the new inverted 
structure could exhibit an even bigger band gap which for the case that the 
initial structure is infiltrated by sihcon gaps, could be of the order of 27%. 

2.4 Applications 

One of the most immediate applications of photonic crystals was to be in op- 
toelectronic devices where unwanted spontaneous emission affected their per- 
formance. On the left of Fig. ^| we show the electron dispersion relation for 
a direct gap semiconductor. Assuming that is embedded in a photonic crystal 
(whose photonic dispersion relation is shown on the right of the same graph) 
it is easy to see that if the PBG overlaps the electronic band edge then the 
electron-hole recombination rate could be inhibited. That is simply because 
the potentially emitted photon will have no place to go! In a semiconductor 
laser this would lead to near unity efficiency into the lasing mode. This theo- 
retical idea was recently implemented "27 in the 2-D band edge micro-laser in 
which lasing occurs preferentially at the 2-D photonic band edge even though 
the emission from the active region has a broad frequency distribution (Fig. I14|l . 

Along the same lines, light extraction from LEDs would be more efficient. 
The main problem there is that only the light emitted in a narrow angle ^ 
manages to escape and the rest is usually trapped within the film and is ei- 
ther absorbed or emerges from the edge of the device. There are several ways 
to partially overcome this, for example using photon recycling, rough surface, 
e.t.c, which could increase the efficiency to 30% but none of them actually alters 
directly the spontaneous emission properties of the device. If a 2-D photonic 
crystal pattern in the form of a triangular lattice of air holes is introduced into 
the semiconductor I^Sj, extraction efficiency of the order of 100% is expected. 
This is simply due to the large band gap that will appear for in plane prop- 
agation forcing the photons to come out of the slab in the vertical directions. 
In addition, unlike the planar microcavity, the efficiency is enhanced over a 
wide range of frequencies and as no resonance or photon recycling is needed, 
the photon lifetime is shorter. The latter of course leads to a reduction of the 
absorption losses and an increase of the modulation speed of the LED. 

We proceed by describing the potential applications of the presence of defect 
modes in a photonic crystal. It was shown that by adding or removing a piece of 

^The efficiency being l/(2n^) where n is the index of refraction for the emitting material. 
For GaAs, n=3.5 and the efficiency is approximately 5 %. 
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Figure 12: (Top): The spiral crystal made of square spiral posts in a tetragonal 
lattice proposed by Toader and John. It is believed that structures like these 
might provide wider and more robust to disorder band gaps |26| . (Bottom): The 
Si crystal obtained by inverting this template. A gap of 23.6% occurs between 
fourth and fifth bands of the photon dispersion relation. The corresponding 
total DOS is shown in the inset (arbitrary units). 



15 



Conduction 



Electronic Baiid ^ 
Gap ^ 



1000 



Electronic dispersion 




I Photonic 
Band Gap 



Valence 



Pliotonic dispersion 



Figure 13: Diagram showing the dispersion relation for both the photon and 
electronic modes in a semiconductor integrated with a PBG material. The EM 
dispersion relation for the PBG structure is shown on the right whereas on the 
left is the electron wave dispersion relation typical of a direct gap semiconduc- 
tor. If an electron for the conduction band were to recombine with a hole from 
the valence band and the PBG overlaps the electronic band edge, the result- 
ing photon would not be emitted as there are no modes available there! This 
will result in the inhibition of the electron-hole recombination rate, seriously 
improving the performance of devices like semiconductor lasers as Yablovonitch 
has first suggested pS) . 
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dielectric material the periodicity of the lattice is locally altered and this can lead 
to the appearance of highly localized modes of light in an otherwise mode free 
surrounding. In particular, removing a small amount of high index material from 
one unit cell (air defect), leads to the occurrence of a localized mode just above 
the top of the lower band in analogy to the acceptor modes in semiconductors. 
On the other hand adding a small amount of high index material to a single 
unit cell (dielectric defect) forces a single localized mode to split off from the 
upper band edge as in semiconductor donor modes. The former case- an air 
defect-is basically a high-Q microcavity with tunable frequency ^. All types of 
applications that involve high-Q optical microcavities could thus implemented. 
The main advantage over usual high-Q superconducting metallic cavities is the 
operation on higher frequencies without almost any losses^. Various devices such 
as frequency filters, atom masers, zero-threshold lasers should be more efficient 
(see Fig. 114b V Also experiments involving single atom-photon interaction in 
the strong coupling regime is important should be easier realizable. 

In a similar way to the introduction of points defects, someone can also 
create line defects in an otherwise perfect structure. Such defects could be used 
as lossless waveguides [21 1221 1^ Hil HE] where light could be guided around 
sharp corners with no reflection or scattering losses (see Fig. I15|l . Also highly 
efficient fibers confining light within a hollow core (a large air hole) in a silica-air 
photonic crystal fiber has been demonstrated 33 . 

Last but not least, we mention the idea of creating tunable PBGs by the 
infiltration of inverted opal or spiral structures with some low index of refraction 
liquid crystal. The potential of controlling the direction of the nematic liquid 
crystal molecules through an external electric field, could alter the optical prop- 
erties of the whole structure and simply provide external control over the width 
of a full 3-D gap |37] . The resulting tunability of spontaneous emission, waveg- 
uiding effects, and light localization may considerably enhance the technological 
value of a composite liquid crystal PBG material over and above that of either a 
bulk liquid crystal or a conventional PBG material by itself. A tunable optical 
microchip which routes light from a set of optical fibers and could be used as 
part of greater optical circuit is shown in Fig. 1161 

We conclude this part by noting that the intense interest in the field of 
photonic crystals in recent years has not solely been due the above described 
"optical circuit" type of applications, where light was treated as a classical 
field. A lot of research has also been towards the understanding of the light- 
matter interactions within a photonic crystal at the quantum level. Several 
novel phenomena, such as, for example, the inhibition of spontaneous emission 
and the creation of "atom-photon" bound states pni^B^HlB51II51IT7lB5l 

l3aisniiniElE3IMlE3I^EZIEaEiEniinilEaE3EJ. and 

also the existence of transparency near a photonic band edge ,69, [73 [73] 

were found in these studies. Describing the effect of inhibition of spontaneous 

*The frequency of the defect mode is an increasing function of the volume of the defect 
area [3 and its energy is exponentially localized (usually within a few lattice constants). 

^It was experimentally shown 1301 that the quality factor increases exponentially with the 
size of the crystal. 
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emission and the creation of "atom-photon" bound states will be the task of the 
section to follow. 

3 Spontaneous emission 
3.1 The influence of the reservoir 

It was not until the 1980's 74 before it was established that spontaneous emis- 
sion is not an intrinsic property of matter over which we have no control on |75| 
but a process greatly dependent on the nature of the surrounding environment. 
We should note here that it was actually 1946 when an effectively overlooked 
paper by Purcell ,76; already had suggested that the spontaneous emission rate 
of radiating dipoles can be tailored by using a cavity to modify the dipolc-field 
coupling and the density of the available photon modes. If the modal density in 
the vicinity of the frequency of interest is greater than that of free space, then 
spontaneous emission will be enhanced, if it is less, it will be inhibited This is 
the Purcell effect. Technological developments, in the form of high quality and 
finesse cavities, extended this to what we know today as cavity QED 77^. One 
of the most characteristic concepts that was introduced at the time was the 
vacuum Rabi oscillations. To illustrate this, assume you have an atom, excited 
to some Rydberg state, placed inside a microwave cavity of high quality factor 
of the order of 10^" and mode spacing larger than the mode width. Assuming 
that the frequency of the transition from the initial Rydberg state to the near- 
est one matches one of the cavity modes, then spontaneous emission is greatly 
enhanced, whereas it is severely inhibited when the atomic transition is detuned 
from all cavity modes by an amount larger than the mode width. In the case of 
exact resonance, the spontaneously emitted photon is reabsorbed by the atom 
then reemitted again and so on, leading to an oscillatory exchange of energy 
between the atom and the cavity radiation field. The operation of the micro- 
maser is based on this effect where a constant exchange of coherent microwave 
photons between the mirrors of a microwave cavity and Rydberg atoms flying 
through |75] . 

From the implementation point of view, the technological expertise required 
to built high-Q cavities has only recently moved from the microwave to the 
optical or near optical regime allowing for the possibility of similar research in 
the optical range 79 . The great obstacle in localizing light in cavities is related 
with the shape of cavity modes. Basically, there are always tails of the usual 
Lorenzian distribution of the cavity lineshape that extend to infinity allowing 
an excited atom located in the cavity to eventually decay to the ground state. 

Extensive research in dielectric structures of the form described in the previ- 
ous section showed that it was possible to fabricate situations where the modal 
density of the electromagnetic field exhibited gaps for a range of frequencies. 
As a qualitative approximation we would say that an excited atom tuned with 
the gap should never decay as there are no modes for its photon to exist. Ini- 
tially lead by John and collaborators and then quickly spreading to many groups 
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Figure 14: Designs for 2-D photonic crystal micro-lasers, (a) The Band Edge 
micro-laser. Stimulated emission (arising from electron-hole recombination) 
from the multiple quantum well active region occurs preferentially at the band 
edge. It has been proved that strong feedback and memory effects can arise in 
that case |2H|- For a real 3-D case this would lead to lasing without a conven- 
tional cavity [391 flU] . The Noda group has realized a precursor to this where 
broad frequency emission occurs preferentially at the band edge f27j (courtesy of 
S. Noda, Kyoto University), (b) Defect Mode micro-laser where a defect mode 
with a localized state of light within the 2-D PBG structure is engineered by al- 
lowing for a missing pore in the 2-D photonic crystal. Stimulated emission from 
the multiple quantum well active region falls mainly into the localized mode 
(courtesy of Axel Scherer, California Institute of Technology). 



19 



S / O is « 9 O * O 

,0 O 0" ^0 o 'T" o 

o o « e ' >5 

c * '* ■ 




Figure 15: Waveguiding in photonic crystals, (a) Right angle waveguide channel 
in a 2D photonic crystal. The propagation of the field modes is through a line 
a defect with no reflection or scattering losses (courtesy of J.D.Joannopoulos 
group, Massachusetts Institute of Technology), (b, c) Other line defect mi- 
crostructures fabricated in macroporous silicon with a lattice constant of 1.5 
mm. The splitting of the modes in (b) could have application in building optical 
devices such as interferometers. In (c), the air holes within the line defect oper- 
ate as reflectors or mirrors thus creating a resonator cavity within the waveguide 
(courtesy of Max-Planck Institute for Microstructure Physics, Halle, Germany). 



throughout the word, an enormous amount of work has been produced dealing 
with questions usually asked in quantum optics, now extrapolated to PBGs, at 
least at their theoretical properties under somewhat idealized situations using 
generalized density of modes for the PBG material [HI ESI EH E3 EEl EH 
EHlESEHIEIlEaESlElESlEniEZlEHlEniEaiMIEHllMlEZlor under more 
involved situations where the local density of modes is used for the description 
of the PBG material [^01 1^ ESI EH ■ The purpose of the sections to follow is 
to review some of the tools used in the above studies and more specifically in 
describing the dynamics of spontaneous emission. The first part describes the 
case where the atoms interact with the free space vacuum and the second with 
structured continuum of a photonic crystal. 

3.2 Two-level atom in free space: Weisskopf-Wigner the- 
ory 

Suppose now we initially prepare our two-level atom in the upper state of the 
doublet and allow it to interact with the modes of the vacuum. The correspond- 
ing Hamiltonian, in the dipole and rotating wave approximations, is given by 

IHOIEI] 

H = hwi\l){l\ + hLU2\2){2\ + hJ^LOkPlPk + n^(|2)(l|/3fc52i(^fc) + h.c.) , (16) 

k k 
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Figure 16: A futuristic depiction of an electro-actively tunable PBG routing 
device. The PBG structure has been infiltrated with an optically anisotropic 
material such as a liquid crystal (indicated by the yellow arrows) which responds 
strongly to external electric fields. Applying a voltage, could alter the optical 
properties of the whole structure opening or closing of the corresponding pho- 
tonic band gap could be achieved. This could lead to the possibility of routing 
light from an optical fiber into one of several output fibers (courtesy of S. John's 
group, University of Toronto). 
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where 



and 

A21 = y 02(x)ex(/)i(x)d^x, (18) 

is the dipole moment for the transition |1) (lower state with wavefunction (f)i(x.)) 
to |2) (upper state with wavefunction 02 (x)). Here, Pk is the annihilation op- 
erator for the field mode k, V is the quantization volume and eo is the electric 
permittivity of free space. Also, fkVj, j = 1,2 is the energy of the j-th atomic 
level. Finally, hujk and denotes the energy and the polarization vector of the 
fc-th reservoir mode. 

Transforming the system in the interaction representation j8U[l81j we obtain 
the interaction Hamiltonian as 

if(^) - hJ2 |2)(l|/3fe52i(^fe)e-''("'=""^^)* + h.c. , (19) 

k 

where huji2 = h{u!2 — i-^i)- Assuming that initially the field modes are in the 
vacuum state |0), the wavefunction of the system can be written in terms of the 
state vectors as 

m)) = a2{t)\2,0) +Y,bk{t)\l,k) , (20) 

k 

where the coefhcients bk represent probability amplitudes of emitting one photon 
of energy Wfe belonging to the vacuum modes |1, A:) (with the atom in state |1)). 
Also 02 (t) represents the probability of the atom being in the excited state and 
the field has no photons. The time evolution of the system is described by the 
Shrodinger equation 

^n^^\m) ^ H^'^im) , (21) 

which using Eq. (|19|l leads to the following set of first order, coupled linear 
differential equations 

ia2{t) = ^.g2i(wfc)e-'''=*6fc(i), (22) 

k 

Mt) = 52*iK)e"^'*a2(i), (23) 

where 5k = ^k — ^12 is the detuning of the emitted photon frequency from the 
atomic transition. We formally integrate Eq. H23|l and obtain 

hu = ~i9*2i{^k) f a2{t')e''''' dt' + bk{0) . (24) 
Jo 
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Taking into account that bk{0) = 0, we substitute Eq. H24|) into Eq. H22I) and 
obtain 

a2{t) = - f dt'a^it') |52iK)pe-'^'=(*-*') , (25) 
•^0 k 

which is still an exact equation. In order to proceed we have to calculate the 
kernel 

^(i"0-El32iK)l'e-^'^'*'*'^- (26) 

k 

The sum over modes is generally, converted to an integral, by including the 
density of modes (states) p(wfc), 

^|g2i(^fc)Pe-^'=(*-*') = ^^ j dnJ^^\gM^,rpMe-'''^'~''\ (27) 

where d^l is the solid angle and a the light polarization. In free space, the 
product \g2i{oJk)\^ p{oj) ~ w is a smooth varying function of uj. Taking this into 
account we can approximate 



dn / |52i(^fe)Pp(c^)e-'("^-"^^'(*"*') ~ (28) 
Jo 

/"OO 

92i{co2i)\^p{i02i) dn e-^(-.-2i)(t-*') ^ (29) 



which eventually gives [SUl |M] 

K{t~t')^:^5{t-t'), (30) 

where 712 = wi2'^/i2i^/(37reo?ic'^) is the famous free space decay rate. This 
approach is the much used the Weisskopf-Wigner approximation. 

The fact that the so called response function K{t — t') is a delta- function 
means that the free space acts as an immediate response reservoir. In other 
words spontaneous emission is dealt as a Markovian process and the evolution 
of the system depends only on the present and not on any previous state of the 
reservoir. Substituting Eq. into (^3)) our initial set of equations become 

d2it) = --^a2{t), (31) 
^bk{t) = 9l^{u;u)e'''''a2{t). (32) 
The above equations can easily be solved with respect to time which gives 

|a2(i)P =e-^-*. (33) 
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3.3 Two-level atom in a photonic crystal 

Assume now that our two-level atom is coupled to the radiation field of a mod- 
ified reservoir. We will specifically assume that the atom is embedded in a 
three-dimensional photonic crystal where the photon dispersion is found to be 
an isotropic one and satisfies the following transcendental equation^*' 



c 

ojfc = - — arccos 
Ana 



Ancos{kL) + (1 - n)^ 
(l + n)2 



(34) 



The Hamiltonian and the wavefunction of the system are still given from Eqs. 
ifT^ . lf^ and Eqs. lf^ . (|^ still describe the evolution with the difference now 
that atomic transition is occurring in the vicinity of a band gap in the density of 
allowed photon modes inside the photonic crystal. As we showed in the previous 
section the excited state amplitude 02 is given by 

02(0 = -/ a2{t')K{t-t') dt' . (35) 
Jo 

In order to calculate the kernel in this case we cannot use the Weisskopf-Wigner 
approximation as the density of states changes very rapidly in the vicinity of 
the atomic transition when located near the band edge. In this case, we must 
perform an exact integration in Eq. H2t)|) . To do this, we observe that in the 
vicinity of the gap the dispersion relation Eq. H34|l can be approximated as 

ujk=ujg + A{k - kof , (36) 

where A ~ ^glk\ (this corresponds to a density of states of the form p(u>) ~ 
0(0; — LOg)/ ^oj — ojg shown in Fig. 1171 with Q being the Heaviside step function 
|41l). The latter combined with Eq. H27|l gives for the response function of Eq. 

K{t-t') = ^ , , t>t', (37) 

with /?^/2 = LL'[2'^|/ii2p/(67reofi-c'^) and Sg — ujg — LO12. Eq. (|T7|) demonstrates 
the non-Markovian character of the reservoir. In contrast to the free space case 
we can see from Eq. H37|l that there is a contribution in the current dynamics 
at time t from previous states of the system at time t' following an inverse 
square root dependence. As we will discuss later on, this is due to the partial 
localization of the emitted photon in the vicinity of the atom where it can be 
re-absorbed and thus affect the atom's evolution again after its initial emission. 

We continue by trying to derive the explicit time dependence of the atom's 
evolution. To do that we need to solve the integro-differential equation (|35(l . 
For this, we first Laplace transform it and obtain 

s~a2{s) - a2{Q) ^ k{s)~a2{s) , (38) 



lOWe remind here that n is the refractive index of the scatterer, a is its radius, and 2a+b = L 
is the lattice constant (see Section 2). 
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where d2{s),K{s) are the Laplace transforms of a2{t) and K{t) accordingly. 
Using Eq. ^ we obtain for the excited state amplitude 0,2(3) 



0,2(3) 



5(s-i(5g)l/2 _ (if3)3/2- 



(39) 



To obtain the dependence in the time domain we invert this using the 
Bromwitch formula: 



a2(t) = 



1 

27ri 



e+ioo 



(40) 



where the real number e is chosen such that s = e lies to the right of all the 
singularities (poles and branch points) of the function 02(5). The inverse Laplace 
transform of Eq. yields gS| 



a2 



(t) = 26iXie''^i+**<'* + 62(2^2 + y2)e'^'='+'^<'' 



where 



Xi 
X2 



Y,b,y,[l-^^(3xp)]e''^^+^'^', 



t/4 



iA+ + A^)e' 
(A+e-^(^/6) - A_e*('^/6))e-*(''/''\ 



1 1 
2±2 



4(5„ 



-,1/2 



27^3 



1 1/3 



(j k;j,i,k^ 1,2,3), 



(i = 1,2,3), 



(41) 

(42) 
(43) 
(44) 

(45) 

(46) 

(47) 



and ^(x) is the error function 83 . As we see from Eq. H42|l . if 5g = then 
Pxf = i(3. This means that the value of (3 given above can be considered as a 
resonant frequency splitting, an analog of the vacuum Rabi splitting in cavity 
quantum electrodynamics [S^. For large times, i.e. for large values of (3t, the 
terms of higher order than {(31)^^"^ [53] can be ignored, and Eq. (|41|l reduces to 



a2(t) ^ 2&iXie^^i+''«* + b2(x2 + 2/2)e''"'+*''* + 77V 



3 J, 



JSt 



{(3t)y^ ■ 



(48) 



As we see from Eq. (I41|l the atomic level splits into dressed states caused 
by the atom and its own radiation field located at the frequencies ujdi — tOg — 
plm(xl) and LLjd2 = ^3 ^ (3Im(x\). It is easy to see using (|42ll and (|45|) that 
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Figure 17: The density of states under consideration. wi2 is the atomic fre- 
quency and LOdiTL0d2 the corresponding dressed states frequencies. These states 
emerge due to the "self dressing of the atom by its own locaUzcd radiation 
field. 



x\ = i|a;ip. This translates to the fact the corresponding dressed state at the 
frequency lOg — I3\xi\'^ is the photon-atom bound dressed state with no-decay. A 
photon emitted by the atom in this state will tunnel for a length scale of a few 
lattice constants before being reflected back and re-absorbed by the atom. The 
other dressed state on the contrary is pushed outside the gap, where the density 
of photon states is not zero and eventually will be responsible for the part of 
the initial excitation that will eventually decay (see Fig. I17II . In Fig. 1181 we 
show the atomic population as a function of time for various detunings from 
the band edge frequency. As was expected for atomic transitions well inside 
the gap {5g = —10/3), the atom remains in the excited state for ever. In this 
case the second term in Eq. (j41|l vanishes which means that no true atomic level 
splitting is present in this case and the oscillations in the dynamics are caused by 
the interference with a "quasidressed" state originating from the branch point 
term. As 6g moves from negative to positive, the other component, located 
at ujg — /3|a;ip becomes important and eventually for 5g — 10/3 is the major 
one forcing the atom to complete decay in the usual exponential manner. We 
note here that this photon-atom bound state is present even when the atomic 
frequency W21 is placed outside the gap where the density of states is not equal to 
zero and that is because of the special singular behaviour around of the isotropic 
model there. For "smoother" cases such as smoothed models for the PBG 
density of modes or a Lorenzian density of modes both the corresponding dressed 
states occur at complex frequencies which for long times leads to complete decay 
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Figure 18: Atomic population on the excited state, P{t) = |a2(t)P as function 
of time for various values of the detuning from the band edge 5 = —10/3 (thin 
solid curve), S = —3.5/3 (thin dashed curve), S = —1(3 (thin dotted curve), 
5 = (thick solid curve), (5=1/3 (thick long-dashed curve), 6 = 10/3 (thick 
short-dashed curve). The time is in units of 1/(3. 

of the population to the ground state. 

4 Conclusions 

In this introductory review article wo have discussed phenomena in photonic 
crystals with emphasis on the inhibition of spontaneous emission of atoms in 
such materials. We started by describing the physics behind photonic band 
gap formation by making a parallel with the known phenomenon of electron 
localization in solids and continued by presenting a simple model of a band 
structure calculation. We briefly presented the state of the art methods in 
fabricating photonic band gap materials. Some of the exciting applications 
in the field of optoelectronics were illustrated followed by noting the basically 
theoretical but really exciting possibility of circuits of light in all an optical 
computing device. In the second part we moved to the quantum regime, i.e., 
the interaction of the quantized electromagnetic field inside a photonic crystal 
with small atomic systems. We discussed the way of treating the spontaneous 
emission of a two-level atom in free space using the Weisskopf-Wigner theory. 
We concluded this introductory review by presenting in detail the modification 
in the spontaneous emission of the same two-level atom when embedded in a 
photonic crystal of a specific type. New phenomena arising by the localization 
of light as the complete inhibition of the spontaneous decay, atom-photon bound 



27 



state and population trapping in a two-level atom were derived and discussed. 
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